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ABSTRACT 

We study the clustering of LRG galaxies in the latest spectroscopic SDSS data releases, 
DR6 and DR7, which sample over 1 Gpc'^/h'^ to z=0.47. The 2-point correlation func- 
tion ^(cr, tt) is estimated as a function of perpendicular a and line-of-sight tt (radial) 
directions. We find a significant detection of a peak at r ~ llOMpc/h, which shows as 
a circular ring in the a — n plane. There is also significant evidence for a peak along the 
radial direction whose shape is consistent with its originating from the recombination- 
epoch baryon acoustic oscillations (BAO). A £,{(J,tt) model with no radial BAO peak 
is disfavored at 3.2a, whereas a model with no magnification bias is disfavored at 
2(7. The radial data enable, for the first time, a direct measurement of the Hubble 
parameter H{z) as a function of redshift. This is independent from earlier BAO mea- 
surements which used the spherically averaged (monopole) correlation to constrain 
an integral of H{z). Using the BAO peak position as a standard ruler in the radial 
direction, we find: H{z = 0.24) = 79.69 ± 2.32(±1.29) km/s/Mpc for z=0.15-0.30 and 
H{z = 0.43) = 86.45 ± 3.27(±1.69) km/s/Mpc for z = 0.40 - 0.47. The first error is 
a model independent statistical estimation and the second accounts for systematics 
both in the measurements and in the model. For the full sample, z = 0.15 — 0.47, we 
find H{z = 0.34) = 83.80 ± 2.96(±1.59) km/s/Mpc. 



1 INTRODUCTION focus on the LRG's anisotropic redshift space correlation 

^ , , /TT^^, N , , , , n function ^(cr, tt), where tt is the line-of-sight (LOS) or ra- 

Luminous red ealaxies (LRG s) are selected by color and , ,. j • ii, ^ rm 

5_i a \ ; J (^j^l separation and a is the transverse separation, i here are 

fT< • magnitude to obtain intrinsically red galaxies in SDSS r ■ ^ ^ i i Am,i n 

V\i JO three sources ot anisotropy at scales larger than 40Mpc/h. 

■ " " ' (Eisenstein etal 2001). These galaxies trace a big volume, . , , ,. ,. i. , , ,i i 

^ o _', It IS well known that peculiar motion dist orts the correla- 

around 1 Gpc h , which makes them ideal for studying c ■ i • n HF~- — Tlino'TK tj- • i n 

^ ' I II 1 ° tion mnction anisotropically flKaiserl 119871 ). It is also well 

clustering on large scales (see lHogg et al.M2005l ). Attention , ^, ^ ■ ^, , , , i -n 

\ — " known that assuming the wrong b ackground cosmology will 

has been paid especially to the baryon acoustic peak around i j , ■ ^ ■ ^ J71 — 1 — o r, i- 1r>■^?^^ uri^ i 

, lead to an anisotropic 4 JAlcock fc Paczvnski 19791 ). What 

a scale of 100 Mpc/h, because of its value as a standard ruler. . , ■ j_ ^ ■ j_t ^ i- i i • i 

^ ^. . , /„„„^N , , . , , IS under-appreciated is that gravitational lensmg also m- 

In Jjjisenstem etal (2005) the baryon acoustic peak was de- , , . , „ , ,, , , ,. 

, ' troduces an a nisotropy or its own to the galaxy co rrelation 

tected m the spherically averaged two-pomt correlation func- ^^^^^.^^ ^ Mat slItoJ[2000bl : IHuI et all B . B . In this 

tion (i.e. the monopole) using LRG s from an earlier SDSS .„ , c -li.- ^ ■ j- ^ i- 

^ , , , ,^ , . r , n n , NT-., nm^T^o paper, we will present evidence oi this lensmg distortion m 

release (about halt the size ot the fanal data). Both 2di'GRS ^, 

, riT^V,ri . , , .r -, , , thc LOS dircctiou. 

and SDSS spectroscopic redshitt surveys have been used 

to constrai n cosmological parameters via the galaxy p ower A distinct focus of this paper in the series is on the 
spectrum (|Eisenstein et al.1 12004 I Sanchez et al.l |2006|), in- baryon acoustic oscillation (BAO) feature in cr — tt plane. 
eludin g i nformation from the baryon acoustic featu re (|Hiitsil In principle, it can be used as a standard ruler to measure 
l2006al (bl: iPercival et all l2007l : ISanchez et aIll2009^ . Photo- both the Hubble expansion rate H{z), via its observed red- 
metric LRG catalogs cover a larger volume and larger den- shift span in the radial direction, and the angular diameter 
sities an d have also been used to obtain cosmologica l con- distance Da(z) . via its observed angular size in the trans- 



sities an a nave also oeen usea to optam cosmologica l con- distance Ua\z) . via its observed angular size m the trans- 
straints (|Padmanabhan et al.|[2007l : iBlake et abllioOTl). verse direction (jBlake fc Glazebrookll2003l : ISeo fc EisensteinI 



In this pap er, as in Papers I ("Cabre fc G aztafiagal l2005h . We will implement this idea to measure H{z). We are 
l2009al ) and II ((Cabre fc Gaztafiaga 2009b., ) of this series, we aided in this endeavor by two effects. Redshift distortions 
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in the LOS direction make the ^ negative on intermediate 
scales (~ 50 — 90 Mpc/h), while magnification bias gives ^ 
a positive boost especially on large scales (ft 100 Mpc/h). 
This combination enhanced the contrast of the radial BAO 
peak with respect to the BAO in the perpendicular direc- 
tion. Because the noise is shot-noise dominated any increase 
in the signal, such as bias or magnification bias, increase the 
signal-to-noise and helps the BAO detection. We will show 
that the shape of the correlation function in the radial direc- 
tion is in good agreement with the predictions. We will also 
validate the interpretation of the feature as baryon acoustic 
in nature by studying its appearance in directions away from 
the radial direction, and performing a parametric fit to the 
monopole. 

It is useful to point out some related earlier work. 
Redshift distortions in the LRG's and quasars at z ~ 
0.55 have been studied using the 2dF-SDSS LRG and 



QSO Survey (2SLAQlCannon fc eta]ll2007l : IRoss et al.ll2007l: 
Ida Angela et al.l l2008l : IWake et al.l boOsE . lOkumura et all 
1 20081 ') measured (.{ayTv) away from the LOS direction in 
the SDSS using a sample similar to ours: they used about 
47000 LRG's over a redshift range of 0.16 - 0.47 while we 
use 75000 LRG's over z of 0.15 - 0.47. None of the earlier 
papers attempted a direct measurement of H{z) from the 
baryon feature in the LOS direction. 

This paper is organized as follows. Section [2] gives a 
summary of the theory, including a brief introduction to 
magnification bias. In section |3] we perform an analysis of 
the clustering signal around the position of the BAO peak, 
both in the monopole and in the a—n plane, especially in the 
radial direction, from which we obtain the Hubble parame- 
ter H(z) in section [l] We end with a conclusion in section 
[S] where we deduce implications for the dark energy equa- 
tion of state, and where we emphasize the need for further 
theoretical work on the nonlinear coupling between peculiar 
motion and magnification bias. 
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2 THEORY 

The 2-point correlation function, (,{r), is defined by the joint 
probability that two galaxies are found in t wo volume ele - 
ments dVi and dV2 placed at separation r (see |Peebleslll98(j ): 

dPi2^n''[l + ^{r)]dVidV2 (1) 

where n is the mean number density of galaxies. 

We can split the distance r into its component along 
the line-of-sight (LOS) tt and perpendicular to the LOS a, 
where = tt'^ + a'^. Azimuthal symmetry implies ^ is in 
general a function of tt and a alone. 

The correlation ^{a, tt) is related to the power spectrum 
by a Fourier transform: 

^(-.-)-/^»W^-"''(§3 (2) 

Note that the use of tt to denote the LOS separation is 
conventional. Hopefully the reader will not confuse it with 
the numerical constant tt as in {2tt)^. 



Figure 1. Top panel: Theoretical ^(cr, tt) with linear redshift 
space distortions convolved with a dispersion model , assuming a 
cosmology and bias (linear and non-linear) as observed in LRG in 
paper I and II. Note the ring around ~ 100 Mpc/h that becomes 
narrow and less prominent in the radial direction. Bottom panel: 
same with magnification bias added (slope = 2 to see clearly the 
effect). The main effect is the boosting of the BAO peak in the n 
direction. In both cases we take into account the bin smoothing 
at 5Mpc/h. 

2.1 Redshift Distortions by Peculiar Motion 

In the large-scale linear regime, and in the plane-parallel 
approximation (where galaxies are taken to be sufficiently 
far away from the observer that the displacements induced 
by peculiar velocities are effectively parallel), the distortion 
caused by coherent infall velocities take s a particularly sim- 
ple form in Fourier space l|Kaiseijll987l ): 

Ps{k) = {l + PdfP99ik)- (3) 

where Pgg{k) is the power spectrum of galaxy density fluc- 
tuation ^g, ^ is the cosine of the angle between k and the 
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Figure 2. Same as Fig[T] with different panels corresponding to different cosmological parameters as labeled in the figures; left panels 
show change with the baryon density Qi,, middle panels with the scalar spectral index ns and right panels with the matter density Qm- 



line-of-sight, the subscript s indicates redshift space, and (3 
is the growth rate of growing modes in linear theory, the 
dimensionless quantity which solves the linearized continu- 
ity equation V.w + (a' /a)/35g — 0, where the prime denotes 
derivative with respect to conformal time. Assuming that 
the galaxy over-density 5 is linearly biased by a factor b rel- 
ative to the underlying matter density 5m, i-e. 5g = bSm, 
and that the velocities are unbiased, the value of /3 can be 
approximated by 



f3' 



(4) 



rsee lHamiltonlll992l . for a review) . After integration in EqO 
these linear distortions in Ps{k) produce a distinctively 
anisotropic ^(cr, tt). Redshift distortions in the linear regime 
produce a lower amplitude and sharper baryon acoustic peak 
in the LOS than in the perpendicular direction because of 
the coherent infall into large scale overdensities. This is il- 
lustrated in top panel of Fig[l] A characteristic feature of 
this effect is a valley of negative correlations (in blue) on 
scales between tt = 50 — 90 Mpc/h, which as we will show is 
in good agreement with our measurements from real data. 
Such a valley is absent without redshift distortions. 

Redshift space distortions on smaller scales (commonly 
called the finger-of-god effect) are often approximated by a 
dispersion model where £, is convolved with a distribution 
of pairwise velocities, parametrized b y a single parameter 
ffy. De tails are presented in Paper II l|Cabre fc Gaztafi agal 
l2009bl ). Note that in Fig.[T] as in all ^{a, tt) contour diagrams 



of this paper, bins of 5 Mpc/h x 5 Mpc/h are used. The 
finger-of-god effect is much suppressed by using bins of this 
size. Readers interested in an analysis that brings out this 
effect rather than suppressing it are referred to Paper II 
(Cabre fc Gaztafi aga 2009b ) . 

FigH] shows how the shape of ^(7r,o-) changes with cos- 
mological parameters: baryon density Qb, initial spectra in- 
dex Us and cold dark matter density flm- Note how the BAO 
ring becomes more or less asymmetric and how this is cor- 
related with the valley of negative correlations (in blue). 

2.2 Magnification bias 

Gravitational lensing inevitably modulates the observed 
spatial distribution of galaxies. To lowest order, there are 
two effects at work. Imagine a number of galaxies lo- 
cated behind some large mass concentration. Dim galaxies 
that otherwise would not have been detected are brought 
into one's sample by the lensing magnification. This in- 
creases the observed number density of galaxies. On the 
other hand, magnification also increases the apparent area, 
which leads to a drop in the observed number density 
of galaxies. The net lensing effect, known as magnifi- 
cation bias, is controlled by the slope of the n umber 
counts (Turner ct al.' 'l984'; 'Webster et al.' 'l988'; 'Fui 
T988; Narayan 1989; Sc hneider. 198a : .Broadhurst et al..,199. 
iMoessner et al.lll99^ 1: 



dlog^oN{< m) 
dm 



(5) 
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function is therefore 
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Figure 3. LRG number count slope s as a function of limiting 
apparent magnitude rur using all SDSS DR6 photometric catalog. 



where N{< m) refers to the number of galaxies in the survey 
with apparent magnitude brighter than m. Note that to 
estimate s we need to know N{< m) for magnitudes that 
are fainter that our spectroscopic hmit of m = 19.2. We can 
use the parent photometric DR6 sample, which is deeper 
but does not have redshifts. Fig. |3]shows the estimation for 
s in the SDSS LRG DR6 galaxies extracted from the full 
(parent) photometric catalog. Because the parent catalog 
is deeper, it also goes further in redshift. Thus a fraction 
of the fainter galaxies in the photometric catalog will be at 
different redshift to the galaxies in our spectroscopic redshift 
samples. Also note that different selection effects go into 
the photometric and spectroscopic samples. Thus the true 
number count slope in the spectroscopic sample is somewhat 
uncertain and might be different for different redshift slicefl. 
Our magnitude cut is around 19.2, where s seems to lie close 
to 1.5, with a potentially large uncertainty if evolutionary 
and/or selection effects turn out to be important. Given that 
variations of s in Fig|3]are of order As ~ 0.5 we consider a 
range around s — 1.0 — 2.0 for our modeling of magnification 
bias. 

The observed galaxy overdensity Sobs is related to the 
intrinsic (or true) galaxy overdensity Sg by: 

SSobs = Sg + Sfj, (6) 

where S^ — (5s — 2)k, with k being the lensing convergence 
which is simp ly a weighted L OS integral of the mass fluctu- 
ation Sm fsee lHui et al.|[2007l . for details) 

Including linear redshift distortions by peculiar motion 
is straightforward: 

Sobs = Sg + S^ + Sv (7) 

where Sv ~ {l+z)H{z)~^dv^/diT with z being the mean red- 
shift, and being the peculiar velocity in the tt direction. 
Squaring the above expression, the net observed correlation 

Note that the number counts here are limited to some particular 
redshift range, e.g. z=0. 15-0.30, and should not be confused with 
the integrated number counts from redshift zero. In other words, 
the slope of the number count s here should best be thought of 
as being related to the slope of the intrinsic luminosity function. 



^obs — £,gg + £,gv + £,vg + £,vv + £,gfj, + £,fj,g + ^fifj (8) 

The velocity-magnification cross-terms are abse nt by virtue 
of Limber approximat ion and linear theory (see iMatsubaral 
l2000bl : lHui et al.ll2008 l. for details). The first term is the true 
or intrinsic galaxy clustering signal. The next three terms 
account for the Kaiser effect. These four terms together are 
the real space analog of what is shown in Eq.Q. The rest of 
the terms account for magnification bias: the magnification- 
magnification term is very small at our moderate redshifts 
of interest; we need focus only on the galaxy-magnification 
cross terms l|Hui et al.ll2007l . [200^ . They are capable of al- 
tering significantly the observed correlation function on large 
scales. Taylor expanding in vr/x, where tt — \xi — X2I is the 
LOS separation between a pair of galaxies and x is the mean 
sample depth, the galaxy-magn ification cross correlation can 
be written as l|Hui et al.ll2007l ): 

C9m(X1, X2, 02) + 5m9(X1, X2, 02) = (9) 
^i/o'^^m(5s -2)(1 + 2)1X1 -X2| 

J (27r)2 """^ ' 

where Pgm is the galaxy-mass power spectrum. Here, Ho is 
the Hubble constant today, Qm is the matter density, z is 
the mean redshift of the sample in question, and x and G 
with subscripts label the radial distance and angles of a pair 
of galaxies. 

The above expression shows clearly the anisotropic na- 
ture of the lensing corrections to the galaxy correlation func- 
tion. A large LOS separation tt = |xi — X2I is clearly favored 
as it should be, since lensing is more effective if the source 
(background galaxy) and the lens (foreground galaxy) are 
further apart. On the other hand, a small transverse sep- 
aration a ~ xl^i ~ ^2! is helpful for boosting the lens- 
ing correction. This makes sense because gravitational lens- 
ing is strongest when the impact parameter is small. The 
magnification distortion of the correlation function hence is 
strongest in the LOS n direction. Note that the intrinsic 
galaxy clustering strength generally drops with separation, 
even as the lensing correction increases with the LOS sepa- 
ration TT. This means for a sufficiently large LOS separation, 
the lensing correction could dominate. 

It is also important to note that for a small trans- 
verse separation a, the integral of Pgm is dominated by high 
wavenumbers, including ones where both the mass ffuctu- 
ations and the galaxy bias are nonlinear. These nonlinear 
effects can further boost the galaxy-magnification correla- 
tion. Indeed, if a were small enough, one might need to worry 
about strong lensing effects. But since throughout our analy- 
sis, we avoid any pairs with a less than 0.5 Mpc/h (because 
of fiber collision issues), we are safely in the weak lensing 
regime. 

The bottom panel of Fig[l] shows the theoretical expec- 
tation for ^(a, n) with magnification bias taken into account. 
This can be compared against the top panel, which has no 
magnification bias. There is a clear enhancement of the BAO 
in the radial direction {a — 0). The effect in other directions 
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is relatively minor. Here, we have used a small scale galaxy 
bias similar to that in LRG galaxies (see Fig|3]) and a num- 
ber coun t slope of s = 2. The large scale galaxy bias is 2. 
See ^Ijain etal] l|2003l ) and lHui et all (,2007 ) for further 
discussions on these numbers. 



2.3 Beyond Standard Theory 

The above discussion of the theory for redshift and mag- 
nification distortions is more or less standard (though the 
anisotropy induced by lensing is not as widely appreciated). 
There are, however, certain limitations to the standard the- 
ory. For instance, the mapping from real to redshift space 
creates nonlinear effects even on very large scales. A sign of 
this is that caustics can occur in redsh ift space even i f the 
real space over density is rather mo dest ()Hui et al.ll2000l ). As 
emphasized bv lScoccimarrol (|2004l '). depending on the statis- 
tics studied, the Kaiser plus dispersion model for redshift 
distortions can be inaccurate even on linear scales. Simi- 
lar concerns apply to magnification distortions when pecu- 
liar motion are taken into account. The expression in Eq. [8] 
ignores velocity-magnification cross-correlation. This is jus- 
tified in linear theory, but not valid once nonlinear effects 
are taken into account. Moreover, for small values of a (ie 
a < 5.5 Mpc/h in our analysis) the lensing effect at large 
TT separations is dominated by the correlation at r ~ a, 
where the clustering is fully non-linear. We need to account 
for these non-linearities if we want to take the full lensing 
contribution into account. 

We have carried out an investigation of these effects us- 
ing numerical simulations. The details are described in Ap- 
pendixlX] The main conclusions are: 1. Eq. (|9]) seems a good 
approximation to the true magnification bias correction to 
the two point function, even after allowing for nonlinear red- 
shift space effects and their coupling with lensing, provided 
that 2. the fully nonlinear Pgm is used in the integral in Eq. 
© i.e. nonlinear galaxy bias is important. The second point 
is crucial and bears repeating. Suppose one is correlating 
two galaxies that are far apart in the radial direction n but 
close in the transverse direction a. One might naively think 
that since the two galaxies are far apart, only linear galaxy 
bias needs to be considered. This is not true once the lens- 
ing effect is included. A small transverse separation means 
the lensing impact parameter is small. One is then sensitive 
to the galaxy-mass correlation on small scales. Recall that 
the lensing magnification is an integral of mass overdensity 
along the fine of sight. The relevant lensing correction 
comes from correlating the galaxy in the foreground with 
the lensing magnification of the background. The dominant 
contribution comes from the portion of the line of sight in- 
tegral that is in fact close to the foreground galaxy. In other 
words, a theoret ical prediction f or using a linear galaxy 
bias (such as bv lHui et al.ll2007h underestimates the lensing 
effect for small cr's, since the true galaxy bias on small scales 
is expected to rise above the linear val ue. This i s illustrated 
in Fig . 131 from the analysis in Paper II l|Cabre fc Gaztanaga 
l2009bl ). The net effect on is described in Appendix \K[ 
Fig. El 

When modeling data, we can use Eq. ((5} with a non- 
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Figure 4. This is the observed LRG clustering bias as a function 
of scale r, defined by b(r) = \J ^obs ('')/5dm (^) , where ^obs is the 
galaxy correlation function and ^dm is the theoretical nonlinear 
dark matter correlation, both in real space. 

linear galaxy bias taken from Fig. |4l Note that this galaxy 
bias, while inferred from data, is strictly speaking not the 
relevant bias to use - what we need is the nonlinear bias 
relevant for galaxy-mass correlation rather than for galaxy- 
galaxy correlation. In our modeling in l3.8.2l we will therefore 
allow an extra overall normalization factor multiplying Eq. 
([9|, which we refer to as A. We will show that the data fa- 
vors a non-zero magnification bias correction 7^ at the 
2 a level, but our estimate for magnification bias A = \ \& 
1.5-sigma smaller from the actual best fit for A. Given the 
mentioned uncertainties, it is possible the real lensing effect 
is in fact larger than our prediction. 



3 ANALYSIS AND RESULTS 
3.1 Data sample 

In this work we use t he recent spectroscopic SPSS dat a re- 
leases, DR7 and DR6 l|Adelman-McCarthv fc etalll2008l '). We 
use the s ame sampl es and methodology here as presented in 
Paper I (jCabre fc Gaztafiag a 2009a ) of this series. LRG's 
are targeted in the photometric catalog, via cuts in the (g- 
r, r-i, r) color-color-magnitude cube. Note that all colors 
are measured using model magnitudes, and all quantities 
are corrected for Galactic extinction following ISchlegel et al.l 
(1998). The galaxy model colors are rotated first to a ba- 
sis that is aligned with the galaxy locus in the (g-r, r-i) 
plane according to: cx = {r — i) — {g — r)/4 — 0.18 and 
C|| = 0.7{g - r) + 1.2[(r - i) - 0.18]. Because the 4000 
Angstrom break moves from the g band to the r band at 
a redshift z ~ 0.4, two separate sets of selection criteria are 
needed to target LRGs below and above that redshift. The 
two cuts together result in about 12 LRG targets per deg^ 
that are not already in the main galaxy sample. The ra- 
dial distribution and magnitude-redshift diagrams for these 
galaxies are shown in Fig.Al and A12 of Paper I. 

We k-correct the r magnitude using the Blanton pro- 
gram 'kcorrect' 0. We need to k-correct the magnitudes 

^ http:/ /cosmo. nyu.edu/blanton/kcorrect/kcorrect_help. html 
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in order to obtain the absolute magnitudes and eliminate 
the brightest and dimmest galaxies. We have seen that 
the previous cuts limit the intrinsic luminosity to a range 
—23.2 < Mr < —21.2, and we only eliminate from the cat- 
alog a small number of galaxies that lay out of the limits. 
Once we have eliminated these extreme galaxies, we still do 
not have a volume limited sample at high redshift. For the 
2-point function analysis we account for this using a random 
catalog with identical selection function but 20 times denser 
(to avoid shot-noise). The same is done in simulations. 

There are about 75, 000 LRG galaxies with spectro- 
scopic redshifts in the range z — 0.15 — 0.47 over 13% of the 
sky. We break the full sample into 3 independent subsamples 
with similar number of galaxies; low z = 0.15 — 0.30, middle 
z = 0.30 - 0.40 and high z = 0.40 - 0.47. The middle sam- 
ple has a lower amplitude of clustering because it includes 
galaxies with a lower luminosity (see Fig.A12 in Paper I). 
Also note that the radial distribution is quite peaked in the 
middle sample (see Fig.Al in Paper I). The resulting clus- 
tering in this sample has a low signal-to-noise and we cannot 
detect the BAO peak in the monopole (see Fig. 31 and A13 
of Paper I). We will therefore concentrate on results from 
the full sample, as well as the low and high redshift slices 
for the BAO detection. 



3.2 2-point correlation 

The two-point correlation function is defined as 



(10) 



where 5{r) = p(r)/p — 1 is the local density fluctuation 
about the mean p = (p) , and the expectation value is taken 
over difi'erent realizations of the model or physical process. 
As mentioned above, the observed correlation function is 
anisotropic. We estimate the correlation as ^(cr, tt), where tt 
is the separation along the line-of-sight (LOS) and a is the 
transverse separation. The absolute separation (in redshift 
space) is ri2 = ^/cr"^ + tt^. 

In practice, the expectation value above is over differ- 
ent spatial regions in our universe, which are assumed to 
be a fair sample of possible realizations (see Peebles 1980). 
A possible complication with this approach is the so-called 
finite volume effect which leads to an estimation bias, some- 
times referred to as th e integral constraint bias (e.g. see 
iHui fc Gaztanagal l| 19991 ) : lBernardeau et all (|2002l )). For our 
samples and scale of interest (ie < 150 Mpc/h) we have 
checked using a large simulation that such a potential esti- 
mation bias is small compared to the errors. We use a sim- 
ulation about 500 times the volume of the SDSS LRG data 
(MICEL7680 with 453 GpcVh^ 2048^ dark matter parti- 
cles and 107 million halos i n a single box; see iFosalba et al.l 
l|2008l ): ICrocce et al.l (|2009l ')') to estimate the "true" correla- 
tion. We then split this large simulation into 216 mocks and 
estimate the mean and error of the 2 and 3-point correlation 
functions in the 216 subsamples (each similar to the SDSS 
LRG samples). We find that the mean agrees well with the 
true value estimated from the full simulation, well within 
the error estimated from the dispersion among the subsam- 
ples. In other words, any possible estimation bias, on our 



scales of interest, is much smaller than the errorbar. The 
ratio defined by the difference between the true correlation 
(from the full MICE7680 simulation) and the mean corre- 
lation in 216 smaller mocks (cut out of the full MICE7680 
simulation) divided by the dispersion (the error) indicates 
how large is the integral constrain bias as compare to the 
error. It turns out to be insignificant, always smaller than 
0.05% of the error in both dark matter and halo mocks. For 
the monopole we also find that the integral constraint bias 
is always smaller than 0.1% of the error for all the scales 
in our analysis, ie s < 150 Mpc/h (see also paper I and III 
where we show several examples of this). 

3.3 Estimation of ^(a, tt) and diagonal errors 

To estimate the correla tion ^(ct, vr), we use the estimator of 
iLandv fc Szalavl l|l993t ). 



DD - 2DR + RR 
'RR 



(11) 



with a random catalog A''^ = 20 times denser than the SDSS 
catalog. We find similar results when we use other similar es- 
timators, such as DD/RR—1 and D D/DR—1 or estimators 
based on pixel density fiuctuations (ie lBarriga fc Gaztanagal 
l|2002t) l. The random catalog has the same redshift (radial) 
distribution as the data, but smoothed with a bin dz — 0.01 
to avoid possible cancellation of intrinsic correlations in the 
data. The random catalog also has the same mask. (We 
will test robustness against variations in the choice of mask 
and the radial selection function below.) We count the pairs 
in bins of separation along the line-of-sight (LOS), tt, and 
across the sky, a. The LOS distance tt is just the differ- 
ence between the radial comoving distances in the pair. The 
transverse distance a is given by Vs^ — vr^, where s is the 
net distance between the pair. We use the small-angle ap- 
proximation, as if we had the catalog at an infinite distance, 
which is accurate until the angle that separates the galaxy 
pair in the sky is l arger than about 10 degr ees for ^(cr, n) 



pair m tne sky is l arger tnan about iU degr ees tor l;(a,Tv) 
(see ISzapudU |2004l ) and iMatsubaral (|2000al ll. This condi- 
tion corresponds to transverse scales larger than cr = 80 
Mpc/h (cr = 165 Mpc/h) for galaxies at z=0.15 (z=0.34). 

There are two sources of error or variance in the esti- 
mation of the two-point correlation: a) shot-noise variance 
which is inversely proportional to the number of pairs in 
each separation bin b) sampling variance which scales with 
the amplitude squared of the correlation. It is easy to check 
that for the size and number density of our sample, the shot- 
noise term dominates over the sampling variance error. DM 
particles have a much higher density than LRG galaxies and 
for them shot-noise is negligible. We can then dilute the DM 
particles to check how and when shot-noise dominates over 
sampling variance. We confirm this using the same simula- 
tion MICE7680 mentioned above. We create out of the large 
box 216 independent mock catalogs with same density than 
the observed LRG galaxies. Fig. [5] shows the error (square 
root of the variance) of the line-of-sight ^ for dark matter 
(black solid line) and halos (red solid line), computed from 
the dispersion between the 216 mocks - this is by definition 
the true error. Note that we dilute DM particles to match 
the observed LRG density (and N(z) distribution) so now 
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Figure 5. In this figure we show difli'erent error estimates for the 
2-point correlation in the line-of-sight direction tt, averaged over 
(T = — 5 Mpc/h. The black solid line and red solid line shows 
the dispersion for dark matter and halos respectively computed 
from 216 mock catalogs (the dark matter mocks are at z=0.3, and 
the halo mocks are at z=0 chosen with a large scale bias of 1.9). 
We also calculate the Jack-knife (JK) error for each mock and 
we plot its dispersion as a shaded region (gray for dark matter 
and orange for halos). For comparison, we also show the JK error 
for the real sample of LRGs (blue dotted line). The dashed lines 
show the error estimates from our analytic model for dark matter 
(black), halos (red) and LRGs (blue) respectively (see text). 

dark matter mocks are also dominated by shot-noise error 
due to the low density. We compare this error against two 
approximations: Jack-knife error (Gray and orange shaded 
regions for dark matter and halos respectively), and an an- 
alytic error-estimate (black dashed line and red dashed line 
for dark matter and halos respectively). The analytic es- 
timate takes the simple form = Qnoisc A^poisson, where 
A^Poiason is simply one over root A'' shot noise where is the 
number of pairs at the separation bin of interest. For dark 
matter mocks, Qnoisc ~ 1 works very well, but as discussed in 
Paper I, the halos or groups do not quite follow a Poisson dis- 
tribution and their relevant Qnoise = 1.4. The l a tter i s con- 
sistent with the findings of others ISmith et al.l l|2007l ) that 
massive halos, by virtue of exclusion zones around them, 
do not have Poissonian shot noise. Fig. [5] shows that both 
Jack-knife and analytic error estimates work fairly well. For 
comparison, we also show the Jack-knife error for the LRGs 
in real data (blue dotted line) and the corresponding an- 
alytic estimate (blue dashed line). The agreement between 
them confirms the validity of our analytic error model, which 
we will adopt for the rest of this paper. Note that the LRG 
errors are in-between the DM and halo errors because of the 
small differences in number density and value of anoisc- 

3.4 The covariance matrix 

We also use the galaxy (halo) mocks to estimate the covari- 
ance matrix between bins in the whole ^(Tr, a) plane. Figure 
[6] shows the covariance for 5 Mpc/h bins in the LOS, nor- 
malized to the corresponding variance (ie diagonal elements 



Figure 6. Covariance matrix (relative to the variance) for LOS 
correlation in 5 Mpc/h bins. Even if the covariance is not negligi- 
ble close to the diagonal, results in our analysis are very similar 
when we assume a diagonal matrix. 

are unity). Full details of this estimation are given in Pa- 
per L The figure shows that at scales tt > 20Mpc/h the 
covariance is quite small. It is smaller than 10% in most of 
the plane, except in the bins closer to the diagonal where 
it can be larger but still lower 30%. At the BAO position 
TT ~ 110 Mpc/h the covariance is lower than 10% for all 
bins. We do our analysis using this covariance matrix, but 
we have checked that results are very similar when we as- 
sume a diagonal matrix (eg see Fig |16p . so that in practice 
the covariance can be neglected. 

3.5 Results for ^(cr, tt) 

In FigOwe see the redshift-space correlation function ^(cr, tt) 
for the complete catalog (z=0. 15-0.47), and for two different 
slices in redshift: z = 0.15 - 0.3 and z = 0.40 - 0.47. RecaU 
from the top panel of Fig.[l]that, without magnification bias, 
the conventional expectation is that one should see a less 
prominent BAO peak in the LOS direction tt than in other 
directions. Instead, we see from the data. Fig. [71 that the 
observed BAO peak actually gets more pronounced along 
the LOS direction, in qualitative agreement with the bottom 
panel of Fig. [T] which includes magnification bias. Indeed, 
we see a very nice ring associated with the BAO peak in the 
data. Note also in the data, we see these valleys of negative 
correlations (blue) on scales of tt = 50 — 90 Mpc/h, which 
are in accord with the predictions of Kaiser distortions (Fig. 

In separating a from tt in the data we have assumed the 
plane-parallel approximation. This introduces a distortion of 
the BAO scale in the perpendicular direction a > 100 Mpc/h 
when TT is small. This can be clearly seen in the plots, where 
there is an artificial concentration of the signal at an angle 
of a few tens of degrees away from the tt = axis, which 
produces an X shape in our a — n diagrams, especially at 
large a's. In reality, this signal originates from smaller an- 
gles. This effect is explicitly demonstrated in our simulations 
(see Fig.A14 of Paper LCabre fc Gaztanaga ,2009a ). In the 
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Appendix of Paper I we show how to correct for this ef- 
fect by removing some of the pairs in the calculation. This 
is not the best possible approach, since it throws away in- 
formation, but it sho ws the origin of th is strange X shape 
feature at large cr's. (|Matsubaral l2000al . see also Fig. 9 in). 
Note, however, that neither the LOS clustering (a = 0) nor 
the monopole, studied here, are affected by this artificial 
distortion. 

We use a fiducial flat model with matter density Sim ~ 
0.25 and Hubble constant h — 0.72 to convert the observed 
redshifts and angles into distances. Unless stated otherwise 
we also assume baryon density ili, = 0.045, and spectral 
index ris ~ 0.98 as the concordance LCDM model {w — —1 
for the dark energy equation of state). 

3.6 The Monopole 

As a first step in studying the BAO, we look at the 
monopole, which is the average of ^ {a, tv) over orientations: 

Ur) = l'jia,n)^ (12) 

where r = \/a^ + TT^ and fi = n/r. Fig.[8]shows the observed 
monopole (solid line with gray area denoting errorbar) in 
different redshift slices. The BAO peak is clearly visible. 
One way to verify its significance is to perform a parametric 
fit and see if a non-zero baryon density fib is required by 
data. We use 4 parameters: Qmh^ , flth^ , and an over- 
all amplitude Amp. Our model is essentially linear theory, 
but includes non-linear effects according to re-normalized 
pe rturbation theory (RPT) at t he BAO peak, as described 
in lCrocce fc Scoccimarrj (|2008l '). Magnification bias is neg- 
ligible for the monopole. As is demonstrated in Fig. 10 of 
Paper I, the ratio of the monopole to the real space corre- 
lation function is constant on scales larger than about 10 
Mpc/h, consistent with the Kaiser model. This is why we fit 
the data with a single amplitude Amp which is supposed to 
account for as, galaxy bias and the redshift distortion boost, 
all rolled into one. To be conservative, only data on scales 
larger than 20 Mpc/h are used in our fit. Covariance be- 
tween different scales is taken into account, using the error 
model developed and tested in Paper I. Magnification bias 



is not included in our fits because its effe ct on the monopo le 
is quite small at our moderate redshifts ()Hui et al.l ()2007h ). 

The resulting constraints are shown in Fig. |9l and the 
corresponding best-fit monopole is shown as a dashed line in 
Fig. IHl The values quoted in the caption are obtained by 
applying singular value decomposition to the covariance ma- 
trix, keeping modes with the smallest 7 eigenvalues, though 
the fits are robust against varying this number. 

In order to calculate the covariance, we have taken 
216 mock catalogs with the same characteristics than LRG 
galaxies. Moreover, we can calculate the JK covariance for 
each mock and compute the mean value over all the mock 
catalogs. This covariance is very similar to the Monte Carlo 
one but smoother so we use the JK one to calculate x^- This 
is extensively explained in Appendix A3 of Paper I. 

There is some tensi on between our best fi t Qb (or Q,bh^) 
and the WMAP5 value iKomatsu et al.l ()2008l ) for the same 
model (flat universe with w=-l). Our best fit values tend 
to be higher. A higher Q.b leads to a more prominent BAO 
peak, as illustrated in Fig. 1101 However, it should be em- 
phasized that in the low redshift slice z — 0.15 — 0.30 and 
in the full sample z = 0.15 — 0.47 the best flt flbh^ is less 
than 2-sigma away from the standard WMAP5 value of 0.22, 
though the discrepancy is larger in the high redshift slice 
z — 0.40 — 0.47. It is also worth noting that there is more 
room for accommodation once more parameters are allowed 
to vary, such as the dark energy equation of state w and the 
ne utrino mass. A more detailed analysis of this is presented 
in ISanchez et aP l|2009h . 

3.7 The BAO Ring 

Recall from Eq. [12] that the monopole receives most of its 
weights from orientations close to the transverse direction 
(i.e. the measure d/i equals sinOdO, where is the angle 
with respect to the radial direction). Let us therefore con- 
sider briefly the reality of the BAO peak in other directions, 
including those close to the LOS or radial direction. In Fig. 
111! we show the signal-to-noise of ^ in the cr — vr plane for the 
redshift slice z — 0.15 — 0.3. This complements the ^(o", tt) 
signal plot in the middle panel of Fig. (TO) The signal-to- 
noise shown in Fig. [Tl]is for each pixel of size 5 Mpc/h by 
5 Mpc/h (the same pixel size is used throughout this paper 
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Figure 8. Measured monopole with errors (solid line with gray area) compared with the best fit model (dashed line), which uses the scales 
20 - 120 Mpc/h (20 log bins). Left panel: corresponds to the full sample (z=0. 15-0.47), where we find = 3.4 (best fit is Q^/i^ = qj2, 
Qbh^ = 0.026). Middle panel: corresponds to the slice z=0.15-0.30, where we find = 1-9 (best fit is Umh'^ = 0.132, U^h'^ = 0.028). 
Right panel: corresponds to the slice z=0.40-0.47 with x^ = 4.8 (best fit is flmh? = 0.124, fl^h^ = 0.04). Allowed values are shown in 
FigEl 





:%^\ fSASi 0.30 025 O.IQ O.I» 0.2C 2^ 10 0.15 0.20 0.35 



Figure 9. Best fit contours (1 and 2-sigma with Idof) for cosmological parameters in a fit to the monopole. Top, middle and bottom 
panels correspond to redshift slices z=0.15-0.47, z=0.15-0.30 and z=0.40-0.47 respectively. The smaller inner (red) contours use a prior 
of flmh'^ = 0.1326 ± 0.0063 from WMAP5. 
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n„=0.25 n.= 0.045 n.= 1.00 n„=0.25 n.= 0.060 n.= 1.30 




Figure 10. A comparison of ^(tr, 7r) in data and models for the z=0.15-0.30 slice. Left panel corresponds to a standard cosmology (low 
Cj, = 0.045 and ns = 1.0) model, while the right panel has a more prominent BAO peak (f2i, = 0.06 and ns = 1.3), which corresponds 
to the best fit monopolc model to the same data. Middle panel shows the data using the same color scheme. 



in all of our ^(cr, n) plots). Note that there is covariance be- 
tween pixels, and so this figure should be interpreted with 
some care (see Paper I). Nonetheless, it demonstrates the 
high quality detection of a BAO ring in the a — n plane. 
The triangle highlights the region n > a, which receives not 
much weight in the monopole, but where the BAO ring still 
shows up nicely. Note that the (S/N)'^ shown is modulated 
by the sign of the signal: the (blue) valley of negative corre- 
lations at TT ~ 50 — 90 Mpc/h - in accord with the predictions 
of the Kaiser effect - are detected with significance as well. 
The overall coherent structure of a negative valley before 
a positive BAO peak (at just the right expected scales) is 
quite striking, and cannot be easily explained away by noise 
or systematic effects. 

Fig. [TT] suggests that there is sufficient information in 
the data to separately constrain the angular diameter dis- 
tance Da{z) and the Hubble expansion rate H{z), as dis- 
cussed in !JT] That is to say, there is in principle the excit- 
ing possibility of determining both quantities, beyond mea- 
suring the combinati on {D\ / H)^^^ from the monopole as 
is customar ily done ( Eisenstein fc etal l2005l : see also Ap- 



pendix A of lHui et al 



2008h . In this paper, we will focus on 



the derivation of H(z) from clustering in the radial direc- 
tion, in part because, as discussed earlier in the plane 
parallel approximation that we have adopted introduces ar- 
tificial distortions to the signal at large a's. This does not 
affect the analysis here which focus on LOS and monopole 
but will have an impact in ^(Tr, cr). Modeling the wide-angle 
effects is possible, and is something we hope to pursue in 
the future. 
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Figure 11. Signal-to-noise ratio in g((T, tt) for the z=0. 15-0.30 
slice. The color scheme denotes (S/N)^ multiplied by the sign of 
the signal i.e. negative values correspond to a negative signal. The 
triangle highlights the region tt > cr, which receives little weight 
in the monopole. 



I), but these variations do not change the location of the 
BAO peak, as will be discussed below. 



On scales larger than 120 Mpc/h we also find some ex- 
cess (i.e. away from null) signal but the amplitude is lower 
and the sign alternates between positive and negative val- 
ues and is less coherent than the positive or negative regions 
at smaller scales. This excess signal on larger scales is also 
found in the monopole for some of the redshift slices. In the 
Appendix of Paper I we look for possible systematic effects 
that could produce this excess. By changing the mask in ex- 
treme ways, it is possible to reduce the amplitude of these 
excess fluctuations on very large scales (see Fig.AS in Paper 



3.8 The Radial Peak 

In this section, we turn our focus to the correlation function 
in the LOS/radial direction. A visual impression of the LOS 
^ in the full sample can be gained from Fig. 1121 [T3l and Table 
[T] which will be explained in more details below. What we 
would like to do first is to check for systematics and test the 
robustness of our measurement. Then, we will examine the 
statistical significance of the BAO feature. 
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Figure 12. This figure shows ^(tt) measured using different ra- 
dial selection functions, masks and datasets, for z = 0.15 — 0.47. 
Details are described in the text. Here, the LOS correlation is 
binned in (tr, tt) pixels of 5 Mpc/h by 5 Mpc/h. Each pixel has 
cr extending from 0.5 Mpc/h to 5.5 Mpc/h (the minimum rr of 
0.5 Mpc/h is imposed to avoid the fiber collision zone). The LOS 
correlation is plotted centered around vr increments of 1 Mpc/h 
- there is therefore an overlap between pixels in the n direc- 
tion. This is purely for presentation purpose, to show how results 
change for different bin center positions. 



3. 8. 1 Tests for Systematics 

There are many checks that need to be made, some of which 
we have aheady mentioned. Here, we provide a complete hst. 

1. Finite volume effects or integral constraint bias, 
namely an estimation bias that results from the galaxy sur- 
vey having a finite size, is examined in 5221 using large nu- 
merical simulations. We find that for the scales of interest, 
the integral constraint bias is totally negligible compared to 
errorbars. 

2. We also employ numerical simulations to test the 
accuracy of our (statistical) error model. This is done in i]3.3l 
(see Fig[S|. We find that our error model agrees well with 
the Jack-knife error and also, in the case of simulated data, 
with the true error estimated from Monte Carlo realizations. 

3. We test for robustness against a different choice of 
the radial selection function. In Fig. 1121 the green zone 
shows the one sigma region of our fiducial measurement from 
DR6. This fiducial measurement, which is used in the rest 
of this paper, uses a radial selection function which is ex- 
actly the observed dN/dz from data but smoothed in bins of 
dz = 0.01. The red solid and dashed lines in the same figure 
show measurements from Data Release 7 (DR7), with dif- 
ferent choices of the radial selection function - the solid line 
uses a radial selection smoothed with dz — 0.05; the dashed 
line uses a radial selection that is not smoothed at all i.e. it 
uses dN/dz straight from the data. We can see that the re- 
sulting ^(tt) is, to within errors, stable against these different 
assumptions about the actual radial selection function. 



4. We test for robustness against a di fferent choice of 
the angular mask. After the release of DR6, ISwanson et al.l 
(2008) provided mask information in a readily usable form, 
translating the original mask fil es extracted fr om the NYU 
Value-Added Galaxy Catalog (iBlanton et al.l 2005). from 
MANGLE into Healpix format (|G6rski et al.ll2005ll . Paper I 
describes how they constructed a survey "mask" for LRGs 
and tested the impact of the mask on clustering measure- 
ments using mock catalogs. Here, we apply the MANGLE 
mask with different completeness in place of our fiducial 
mask to DR6, and the result is shown as the gray zone in 
Fig. 1121 which encompass the whole range of possible com- 
pleteness above zero. Our measurement appears to be stable 
against variation in the mask. 

5. We test for the stability of our measurement when 
we look at different subsets of the data. The triangles in Fig. 
[12] show ^(tt) measured from a set of DR6 LRGs with a dif- 
ferent magnitude cut. The dotted line shows ^(tf) measured 
from the north stripe of DR6. Finally, the red lines show 
measurements from DR7, which is about ~ 17% larger than 
DR6. Again, they are all consistent with each other to within 
errors. 

In Table [1] we list our best estimates and errors for 
the measurement of ^(tt, ct) along the LOS for our default 
DR6 mask. This corresponds to the green shaded region in 
Fig. 1121 Because errors are correlated on scales smaller than 
5 Mpc/h we give the measurements and errors in the Table 
smoothed with a top- hat window of total width 5 Mpc/h, 
this makes the analysis quite insensitive to the choice of the 
central bin position. The resulting covariance between these 
5 Mpc/h bins is negligible in practice (see FigO and the 
different points in the figure can be treated as independent. 
See also Appendix 2.4 in paper I for an extensive treatment 
of the covariance. 



3.8.2 Model fitting: shape constraints 

We next want to check how well the data agree with the 
theoretical model presented in previous sections. In particu- 
lar, we are interested in assessing the statistical significance 
of any BAO detection. This is illustrated in Fig llSI The ob- 
served correlation function (points with errorbars) is shown 
in radial bins of Att = 5 Mpc/h and 0.5 < a < 5.5 Mpc/h. 
We fit the data with a model consisting of the following in- 
gredients: the real space power spectrum is controlled hy Qb 
and Q,rn with an amplitude parametrized by has, where b is 
the linear galaxy bias; redshift space distortion is described 
by a Kaiser model parametrized by /3 (nonlinear pairwise 
velocity dispersion is kept fixed at 400 km/s but its effects 
are important only for scales less than vr < 40 Mpc/h); 
the magnification bias correction is linear in n (Eq. [9]), 
parametrized by an overall normalization A (as described in 
H2.3|l . Note that we ^4 = 1 corresponds to s = 1.5. As we al- 
lowed A to varied this corresponds to degenerate changes 
in s or other normalization factors. The parameters 57 _g, 
f^m, bas and /3 are constrained by the observed monopole 
and quadrupole of ^(a, tt) (Paper I). We therefore impose 
the following (1 — a Gaussian) priors: fis — 0.044 ± 0.003, 
= 0.245 ± 0.020, bas = 1.56 ± 0.09 and /3 = 0.34 ± 0.03 
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Table 1. Our best estimates of the 2-point correlation 4(7r,(T) 
in the LOS for 0.5 < cr < 5.5 ( Mpc/h). In the radial direction 
the correlation is binned with Att = 5 Mpc/h and we take the 
average correlation within that bin. The covariance between bins 
is negligible when using this binning. This is for z = 0.15 — 0.47. 



(the prior on n_g is from WMAP5, while the rest is from 
Paper I). We hold fixed h = 0.72, erg = 0.85 and Us = 0.96, 
which are mostly degenerate with the other parameters, but 
verify that relaxing them in a manner consistent with cur- 
rent data does not alter our conclusions significantly. We al- 
low A to vary anywhere between —1 to 5 (flat prior). Lastly, 
we also allow for a shift Dr in the radial scale vr to ac- 
count for the fact that we do not know the true value of 
E{z) = H{z)/Ho to convert the measured redshift into the 
radial distance tt. Note that we can take this to be indepen- 
dent of Qm in H{z), even for a flat universe, as this would 
correspond to variations in the DE equation of state w which 
does not alter the shape of the correlation on linear scales. 
We define 



Dr 



Hfjz) _ Ef{z) _ ^0.25(1 + 0.75 



(13) 



H{z) E{z) E{z) 

where Ef{z) corresponds to the fiducial fiat LCDM cosmol- 
ogy with Q,n = 0.25 used in our analysis. 

The best fit model (to the full sample) is shown as a 
solid black curve in Fig. 1131 Here, the range of scales on 
which the fit is performed is tt = 40 — 140 Mpc/h. The 
best of the fit is — 8 for 14 degree of freedom (as la- 
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Figure 13. Points with 1-sigma errorbars show the measured 
LRG correlation (,{a, vr) along the LOS (for 0.5 < cr < 5.5 Mpc/h) 
in the full SDSS sample (z=0.15-0.47) as given in Table [T] The 
solid black line is the best-fit model as described in text. The 
dotted black (straight) line shows the contribution from magni- 
fication bias. The dashed black line with no-BAO feature is a 
no-wiggle model. The dashed black line with a BAO feature is 
a no-magnification-bias model. The for the best-fit, and the 
Ax^ for the two dashed lines are given at the top of the figure. 
The dot-dashed blue lines show the 1 — cr range allowed by a shift 
Dr in the radial scale tt. 



beled in the figure), which is reasonable (gives a probability 
of Prob — 87%). We can use this model fitting procedure 
to address the reality of the BAO feature: is it supported 
by data? A model wi th no BAO (generated according to 
lEisenstein fc Hulll998h leads to a Ax^ = 10 away from the 
best-fit i.e. it is ruled out at the 3.2 a level. 

It is also interesting to ask to what extent the magnifi- 
cation bias correction is required. We find the magnification 
bias normalization A — 4.5 ± 2.2. Another way to put it is 
that a model with no magnification bias leads to a Ax^ = 4 
away from the best-fit. The data favor a non-zero magnifica- 
tion bias correction at the 2 a level. Our estimate for magni- 
fication bias, yl = 1, is 1.5-sigma away from the actual best 
fit. Although this is not a very significant deviation, it might 
hint that the real lensing effect is larger that the prediction 
in the Appendix. This could be caused by a slope s > 1.5 or 
by the combination of non-linear redshift space distortions 
and non-linear bias, which are hard to model. 

We should also mention that a simple ^ = model car- 
ries a Ax^ = 6, and so it is only slightly disfavored compared 
to the best-fit model. However, a ^ = model is clearly not a 
good description of the data once we include smaller scales. 

The radial shift parameter Dr is constrained to be: 



Dr = 0.998 ± 0.037 (Icr) 



(14) 



when we marginalize over the other parameters. The best fit 
value of Dr is very close to our fiducial LCDM model with 
f2m ~ 0.25. The ±l-sigma range (i.e. models with Dr — 
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Figure 14. Joint constraints on the radial shift parameter Dr 
and /3 when no prior is assumed on (3. The solid curve shows the 
Ax^ = 1 contour. The vertical lines delimit the 1 — a region on 
P from the observed quadrupole and monopole (Paper I). 

1.035 and Dr = 0.961) are shown in Fig llSI as blue dot- 
dashed lines around the best fit model (black solid line). 
The corresponding value of H{z) is shown in the top entry 
of Table [3] In principle one could find stronger constraints 
in Dr by fitting to smaller n scales, but this will be sensitive 
to the modeling of nonlinear redshift space distortions (i.e. 
finger-of-god). 

Next, we investigate the robustness of our model fits. 
We have systematically explored the constraints on Dr when 
the priors are relaxed. An example is shown in Fig. 1141 Here, 
we remove the prior on l3. The central value for Dr remains 
remarkably robust. The errorbar on Dr does increase, as ex- 
pected, but it is not unduly large. The best-fit /3 from the 
LOS data (with no prior on /3) lies outside the 1 ~ a region 
from the quadrupole and monopole data, but the errorbar 
from LOS data alone is quite large and there is overlap be- 
tween constraints from the two different kinds of data. 

We have also investigated what happens if we fit to a 
different range of scales tv — 70 ~ 140 Mpc/h. The best fit 
model has = 8 for 8 degrees of freedom. A model with no 
BAO has a Ax^ ~ 10 from the best fit, and a model with 
no magnification bias has a Ax^ = 4 for the best fit. The 
results are therefore quite similar to those from tt = 40 — 140 
Mpc/h, but the errors on Dr in Eg 1141 increase by almost a 
factor of 2 for TT = 70-140 Mpc/h, indicating that the H{z) 
constraint is partially driven by the shape of the correlation 
rather than by just the position of the radial BAO peak. 
We will show in next section how to measurement the BAO 
peak position with independence of the shape, which will 
provide an absolute measurement of H{z) (independent of 
h). 

Fig llSI shows the same analysis for the nearby redshift 
slice z = 0.15 — 0.30 which has lower density and much larger 
errorbars. The qualitative results are nevertheless similar: 



Figure 15. Same as Fig[T3]for z = 0.15 - 0.30. 

the model fits the data well, the no-BAO model is ruled 
out with Ax^ ~ 10 and the model without magnification is 
ruled out to Ax ~ 6. The value of Dr is close to unity but 
constraints on Dr increase to about 7% as indicated by the 
dot-dashed lines in the figure. 

For a global view of how well the model and the data 
agree, see Fig. 10 in paper I, where we plot the data plus 
model in the a — n plane for the full sample. Visually, it 
looks quite good. Doing a detailed overall fit to ^(cr, vr) is 
compli cated in practice because of the small scale modeling 
(eg see lCabre fc Gaztaiiaga|[2009lJ ). We leave this task for a 
future analysis and focus here in the line of sight data. 

3.8.3 The BAO Feature 

Next, we would like to focus on the BAO feature in the LOS 
around 110 Mpc/h, and assess its statistical significance 
without recourse to model fitting like in the last section. 
We compute the against a null signal i.e. x^ ~ (S/N)^ 
where S is simply the observed ^ itself. The result is plotted 
in Fig |16l where we show results for the full sample and for 
the low and high redshift slices. Dotted lines use the full 
covariance in Fig[S] while solid lines neglect the covariance 
0. The difference is small reflecting a small covariance in the 
radial distance. The bin width is 5 Mpc/h, so there are 2-3 
independent measurements across the BAO peak. A peak is 
clearly detected except from the intermediate redshift slice, 
z=0.3-0.4, which is not plotted here (x^ just fluctuates inside 
the 2-sigma region). Errors are bigger in this z=0.3-0.4 slice 



Suppose we label the pixels by i. At pixel i, the ^i^" 
glecting covariance is ^f/Ca, and the including covariance 
is CiCi^^Sj • Here, is the signal at pixel i, and dj is the 
error covariance between pixel i and pixel j, and Ca is simply 
the diagonal variance at pixel i. 



© 0000 RAS, MNRAS 000, 000-000 



14 E.Gaztanaga, A.Cabre, L.Hui 




Figure 16. In this figure we plot the significance of the detection (away from zero) as = (S/N)'^ without covariance (soUd Une) 
and including the covariance (dotted line) for bins of 5 Mpc/h (incrementing in steps of 1 Mpc/h). We also plot the results for a more 
conservative mask and a smoother selection function (dashed lines). Each panel corresponds to a different redshift slice, as labeled in the 
figures. The gray zones indicate a l-cr, 2-cr, 3-0" or 4-(7 detection 



because there are fewer number of pairs in the radial direc- 
tion at these scales. This sub-sample also has a lower galaxy 
bieis because it includes less luminous galaxies (see Paper I), 
hence a lower signal and, consequently, less signal-to-noise 
(recall we are shot noise dominated). 

To assess the effects of possible systematics, we have 
repeated the analysis using an angular mask with 10-20% 
less area (and fewer galaxies), which are safely inside "good" 
plates, and with a radial selection function which is smoother 
(dz=0.05) than our default value of dz=0.02 (see Appendix 
in Paper I). The results are shown as dashed lines in Fig |16l 
The BAO peak remains significant in both the low and high 
redshift slices, but less so for the full sample. We should 
emphasize, however, this last statement on the full sample 
only pertains specifically to the (S/N)'^ right around the 
BAO feature at 110 Mpc/h. When we carry out a model 
fitting to the full sample over a wider range in scales (e.g. 
40 — 140 Mpc/h), we find that a no-BAO model is in fact 
always disfavored by more than 3(7, regardless of the mask 
or radial selection function (see i]3.8.2l and Fig llSfl . 



4 A DIRECT DETERMINATION OF H{Z) 

Our discussions above suggest we could attempt to infer 
H{z) from the LOS data in two ways. One is to carry out a 
model fitting as in i l3.8.2l Constraints on the shift parameter 
Dr can be translated directly into constraints on H{z). The 
result we find is shown as the first entry in Table |3l The 
statistical error comes from the marginalization explained 
in the above section, while the systematic error correspond 
to the difference in Dr when we use the different versions of 
the data shown in Fig |12l We call this the "Shape Method" . 
Note that this method gives H{z) in units of -ffo and it re- 
lies on not just the BAO feature around 110 Mpc/h, but 
also the shape of the correlation function on smaller scales. 
Such a method is very analogous to how {D\/ HY^^ is in- 
ferred from the monop ole data in the discovery paper by 
lEisenstein fc etall (|2005l ). 



Sample 
z range (mean) 


''"BAO i ""st ± ^aya 

Mpc/h 


^ZBAO ± Tst ± a ay a 


0.15-0.30 (0.24) 
0.15-0.47 (0.34) 
0.40-0.47 (0.43) 


110.3 ± 2.9 ± 1.8 
110.5 ±3.6 ±2.1 
108.9 ±3.9 ±2.1 


0.0407 ±0.0011 ±0.0007 
0.0428 ± 0.0014 ± 0.0008 
0.0442 ± 0.0015 ± 0.0009 



Table 2. The BAO fiducial scale rgAO i^i the LOS direction cal- 
culated with a flat reference Hj-^f{z) cosmology of £7m = 0.25, 
for three redshift slices: in parenthesis is the respective pair- 
weighted mean redshift, and Oat and cTaya are the statistical and 
systematic errors on r^AO- The direct Azbao measurement is 
shown in the third column, which relates to the fiducial scale as 
Azbao = ^BAO^'^ref (2)/c and is independent of the value chosen 
for H^^f{z). 



4.1 The Peak Method 

The second method to measure H{z) is less model depen- 
dent and focus entirely on the BAO feature around 110 
Mpc/h. We find the peak location, and use that as a stan- 
dard ruler to measure the radial distance. We call this the 
"Peak Method" . We will compare both methods in this sec- 
tion. Note that in the shape method we will only constraint 
H{z) I Ho while in the peak method we will get H{z) directly. 
Operationally, for the peak method, we use the relation: 



H(z)tiuc 



tbao 



H{z)rcl 



(15) 



rwMAP 

where rwMAP = 153.3 ±2.0 Mpc (Table 3 in lKomatsu" et al.l 
(|2008! )) is the comoving acoustic scale inferred the cosmic 
microwave background ( WMAP5) , and tbao is the apparent 
BAO scale inferred from data using the fiducial expansion 
rate H{z)icf ~ //o\/0.25(l + z)'^ + 0.75 to convert redshifts 
to distances. 



4.2 Statistical Error 

The peak measured in the radial direction is broad, just as 
expected from the BAO signal (see Fig |13p . We will take the 
BAO position to correspond to the location of the maximum 
in the (S/N) around the broad BAO peak, and we will use 
Monte Carlo simulations to estimate the associated errorbar. 
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Figure 17. Lower panel: dashed line and dotted lines show the 
mean and dispersion for a model of the correlation function in the 
radial direction for the redshift slice z=0.15-0.30. This model uses 
a smoothed version of the real data as the baseline (the mean). 
The continuous line shows one of the 10* Monte Carlo realiza- 
tions, which have noise added in a way that accounts properly 
for covariance seen in the data. Top panel: = i^/^)^ in the 
model (dashed line) and in the same one Monte Carlo realization 
(continuous line). 



The resulting fiducial BAO scale location tbao are shown 
in Table [21 

The Monte Carlo simulations are generated as follows. 
We model the signal as a smoothed version of the data. We 
use a top-hat smoothing of 5 Mpc/h width and radial bins 
centered in IMpc/h intervals. Such a model is shown as a 
dashed line in Fig |17l for the sample z=0. 15-0.30. Given the 
model, we then add noise with the same covariance and bin- 
ning as in the real data. A detailed discussion of our noise 
model, and the tests we have applied to it, is given in H3.3I 
Note that while we have chosen pixels of 5 Mpc/h in size, 
the data actually contain information on finer resolution. In 
our analysis and in our simulations, the pixels are in fact 
incremented in the n direction by steps finer than 5 Mpc/h 
to facilitate the search for the BAO peak, and the full co- 
variance between these overlapping (and non-overlapping) 
pixels is taken into account. 

Fig llTI shows one realization of such simulations. For 
each realization, we infer the BAO peak location using ex- 
actly the same method as applied to real data. We repeat 
this 10* times and estimate the distribution of the inferred 
BAO peak location. The errorbar is given by the difference 
between the inferred peak location and the input BAO scale 
in the model. The distribution of fractional error in the BAO 
location, and the distribution of peak = (S/N)'^ values, 
are shown in Fig llSl for the same z = 0.15 — 0.30 sample. We 
find an rms fractional error between 2.6 and 3.5%, depend- 
ing on the redshift sample used. The resulting errorbars in 
the peak location are shown as ast, the statistical error, in 
Tabled 



Figure 18. Lower panel: Distribution of relative values in the 
BAO peak location in 10* Monte Carlo simulations, taking a 
smoothed version of the data as the mean for the redshift slice 
z=0. 15-0.30. Top panel: Distribution of peak values of x^- Here, 
= (S/N)^ with S being the difference between ^ and zero i.e. 
this is against null detection. 



4.3 Systematic Error 

There are two sources of systematic errors. One arises from 
the measurement process, and the other is theoretical in 
nature. Let's discus s the theoretical one first. Recently, 
ISanchez et al.l (|2008l ) undertook a thorough investigation of 
systematic effects in the determination of the sound hori- 
zon from the galaxy correlation function. Their Fig. 2 shows 
that at Sim ~ 0.25, the peak in the correlation function sys- 
tematically underestimates the true sound horizon by about 
1.5% (this error increases to 2% for ~ 0.2 and reduces 
to 1% for Q,m = 0.3). One can understand this result by 
modeling the monopole correlation as a sum of a power-law 
and an Gaussian peak center on the BAO position. As we 
change the power-law index and the relative amplitude of 
the peak the maximum of the correlation can shift with re- 
spect to the maximum of the Gaussian. This effect should 
be different (and weaker) for the LOS, because the under- 
laying correlation is flatter (and in fact slightly increasing 
with scale). On the other hand, magniflcation bias tends to 
move the peak b y a similar amount in the opposite direction 
(iHui et al.ll2007l ). Correcting for these biases is in principle 
possible, but requires accurate LOS modeling which is not 
currently feasible. 

A larger source of systematic uncertainty comes from 
the model used in the Monte Carlo simulations in the pre- 
vious section. If instead of the data or a smoothed version 
of the data, we take our best fit theoretical model in section 
13.8.21 as the mean input to the Monte Carlo noise realiza- 
tions we find a much larger error in the BAO position, of 
size 4 — 8% depending on the different assumptions. The 
reason for this is that the model has a lower amplitude than 
the data around the peak and models are quite degenerate 
to uncertainties in the cosmological parameters. Even when 
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Shape Method 


redshift 
range 




km/s/Mpc h72 


0.15-0.47 


0.34 


83.87 ± 3.10 ± 0.84 


Peak Method 


redshift 
range 




H(z) ± Ust ± asys 
km/s/Mpc 


0.15-0.30 
0.15-0.47 
0.40-0.47 


0.24 
0.34 
0.43 


79.69 ± 2.32 ± 1.29 
83.80 ± 2.96 ± 1.59 
86.45 ± 3.27 ± 1.69 



Table 3. The inferred H{z) with its associated errors, ie sta- 
tistical and systematical, for each redshift slice. The top entry 
corresponds to fitting the shape of the LOS correlation to models 
and marginalized over cosmological parameters (note that this is 
in units of Hq = 72^72 Km/s/Mpc). Bottom values use the Peak 
Method explained in Section |4. II which is model independent. 



this is not a statistically significant departure, the smaller 
amplitude in the model results in a large degradation in the 
recovered error in the peak position. We can of course choose 
models that have higher peak amplitude, and that are still 
allowed by the data. Those models will provide artificially 
small errorbars. 

We can account for this systematic deviations in a 
model independent way by noticing that the observations 
must in fact result from the true model plus a realization of 
the noise. Thus, if we add noise to the observations we will 
in fact explore, by definition, all possible models that can 
agree with the data. We can use these new sets of models 
(data-|-noise realizations) as a new baseline to our Monte 
Carlo simulations and explore how much different the error 
can be when we allow for a variation in the baseline. We 
have run 500 such data+noise realizations and have run 10'' 
Monte Carlo simulations for each of them (as in previous sec- 
tion) to find the distribution of possible errors. The result 
is shown in Fig. [TH]for the z=0.15-0.30 slice. The variance 
in the resulting distribution provides us with an estimate of 
the systematic uncertainty in our error determination. 

The statistical ast error is taken to be the variance when 
we use just the data as a based model rather than the mean 
in the distribution of errors from data-|-noise, which is larger 
because the error has been added twice. 

An additional source of our systematic error is esti- 
mated by taking the difference between results obtained us- 
ing different angular masks and selection functions (see Fig. 
[T2]and associated text). Note that our results do not change 
significantly when we use different masks in DR6 or even the 
DR7 data, which became available after the first version of 
this paper. This is then added to the above systematic error 
to produce Oeys given in Table[51 Errors are then propagated 
to H{z) in Table El 

One additional source of systematic error might have 
occurred to the reader. Our estimate of the BAO scale is 
not strictly from the LOS direction: as discussed in i]3.8l we 
have used an array of pixels extending in the tt direction but 
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Figure 19. Distribution of errors in the peak positions using 500 
different baseline models (each model being a realization of the 
data plus noise). For each baseline model we do 10* montecarlo 
simulations to estimate the corresponding distribution and error 
(ie as in lower panel of Fig |17ll . The variance in this distribution 
gives us an estimate of the error in the error. This is for the 
redshift slice z=0.15-0.30. 



with (7 = 0.5 Mpc/h to 5.5 Mpc/h, instead of being centered 
at (T = (to avoid the fiber collision zone). This amounts 
to approximating the BAO scale r = \/tx^ + a'^ ~ tt. The 
lowest order fractional correction is (a/7r)^/2 4 x 10~* 
which is negligible. 



4.4 Results on H{z) 

Table |3] summarizes the different H{z) estimations and 
corresponding errors. The values of H{z) from the "Peak 
Method", ie using Eq.UHl), for the low (z=0.24) and high 
(z=0.43) redshift slices are independent of each other be- 
cause they come from galaxies that are far apart. These 
values are the main result of this paper, as they are both 
model and scale independent (see previous section). In con- 
trast, the values from the "Shape method" have been shown 
to depend on the range of scales used in the fit. The main 
virtue of this "shape method" is to show that the LOS corre- 
lation agrees well with the BAO prediction, which indicates 
that measured peak does in fact correspond to the BAO. 

The values of H{z) from the peak method are in- 
dependent of and in excellent agreement with the pre- 
dictions for H(z) according to the current combined con- 
strain ts from WMAP5, sup ernovae la and the monopole 
BAO iKomatsu" et al.l (|2008r ). assuming a fiat LCDM cos- 
mology (ie w=-l). In Fig |20l we c ompare our estimates with 
Ho = 70±5 km/s/Mpc from HST lvan Leeuwen et all (|2007l ) 
at z ~ 0.03 and with the H{z) inferred from WMAP5 model- 
ing with Ho = 71.9±2.6 km/s/ Mpc and D.m = . 258± 0.030 
for a flat universe LCDM (w=-T) lKo"matsu et al] (120081 ). For 
a flat cosmological constant dominated model, our measure- 
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Figure 20. H(z) obtained using the "Peak Method" over three 
different redshift slices (small dots with errorbars) compared to 
the HST value at 2 ~ 0.03 (big dot with errorbar). We also plot for 
comparison the best value for H(z) from WMAP5 + SNIa + BAO 
monopole (solid line) and its associated allowed region at 2-sigma 
level (dashed lines), assuming a flat w=-l LCDM cosmology. For 
reference, the two dotted lines show the values of H{z) in open 
and closed models without cosmological constant. 



defined to give q = for tlie EdS (fiat universe witliout 
Lambda) Cosmology. Using our two independent values of 
H{z) at z = 0.24 and z — 0.43, we find at intermediate 
redshift, z=0.34: 

q[z ~ 0.34] = 0.93 ±0.33 (18) 

If instead we compare our value for tiie full sample z — 
0.15-0.47 at mean z = 0.34 with tiie HST value Hq = 70±5 
km/s/Mpc at z ~ 0.03 we find: 



g[2~0.18] =0.82 ±0.30 



(19) 



Both results provide independent evidence for cosmic 
acceleration (ie away from EdS model). The mean value 
seems to increase with redshift which is contrary to all ex- 
pectations (cosmic acceleration should turn into deceleration 
at high redshifts) , but this trend is not significant given the 
errorbars. The interest of these results is that for the first 
time we have direct (geometrical) measurements of cosmic 
acceleration at different redshifts. Fig |21l shows the differ- 
ent estimations of q{z) for the 6 possible combinations of 
the 4 values of H{z). Note that these estimates are not in- 
dependent. Current errors are larger than typical evolution 
expected in this redshift ranged for different cosmological 
models (shown as dashed lines) so it is hard to conclude 
more with this data. 
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Figure 21. Estimates for cosmic acceleration q{z) based on the 
different pairs of H{z) measurements shown in Fig |20l Dashed 
color lines show values for flat (top lines) and open models with 
f2,„=0.20 (red), 0.25 (green) and 0.30 (blue). Dotted line corre- 
spond to EdS model (n^ = 1.0 and = 0). 



ments extrapolate to 

Ho = 71.83 ± 1.55 (±1.03) km/s/Mpc 



(16) 



if we use = 0.245 ± 0.020 obtained from fitting ^(o-, tt) 
to the same data in Paper I. 



We can estimate the acceleration parameter 
q[z] = dlni?/dlna±3/2. 



(17) 



5 DISCUSSION & CONCLUSIONS 

In this paper, we have studied in detail a peak in the LRG 
correlation function, in the monopole, in the circular ring 
on the a — 77 plane, and especially in the LOS direction. Its 
location around 110 Mpc/h is consistent with the interpre- 
tation that it originates from baryon acoustic oscillations. 
Its significance can be assessed in several ways. We have 
performed a parametric fit to the monopole: the data re- 
quires a non-zero fij, with high significance. In fact, we find 
that data appears to favor a value that is slightly higher 
than the standard WMAP5 value. However, the difference 
is less than 2 — cr at both the full sample and the low redshift 
slice z = 0.15 — 0.30, but is larger at the high redshift slice 
z = 0.40—0.47 (see Fig. [9l). We have also found this tendency 
to high in paper HI l|Gaztafiaga et al"]|200 8') of this series, 
where we analyze the BAO scale in the 3-point correlation 
function, which tests non-linear growth of perturbations, in 
contrast to the 2-point correlation function, that tests linear 
perturbations. Plots of the signal and the signal-to-noise of 
^ in the a — tv plane (Fig. [10] and lllfl provide a reassuring 
view of the reality of the signal: a negative valley (in blue) 
of negative correlations at tt ~ 50 — 90 Mpc/h together with 
a ring (in red) of positive correlations at ~ 110 Mpc/h. The 
BAO peak in the LOS direction is pronounced, and detected 
with significance (Fig. I13l andll6 p . 

It would be hard to explain such a peak in the 2-point 
monopole, in the plane a — n, in the LOS and in the 3- 
point function with systematic effects. We have checked for 
possible systematics by varying the angular mask and the 
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radial selection function, and find the detection robust us- 
ing either DR6 or DR7 SDSS data (see ^X8A} . The data 
on scales tt = 40 — 140 Mpc/h is well fit by a model that 
includes linear redshift space distortions, magnification bias 
and BAO. In fact, a model with no BAO is disfavored by 
the data at the 3ct level, while a model with no magnifi- 
cation bias is disfavored at the 2a level. Our detection of 
the BAO modulation in the LOS direction is helped by the 
combination of two effects: redshift space distortions make 
^ negative on scales of 50 — 90 Mpc/h while magnification 
bias gives a significant positive boost on larger scales. 

We thus have significant evidence that the LOS corre- 
lation on scales 40 < vr < 140 Mpc/h reproduces well the 
expected BAO signal. Consequently we use the LOS data to 
infer H{z) in two ways which we have labeled the shape and 
the peak method. In the shape method we have used the 
shape of the LOS correlation to fit H{z)/Ho marginalized 
over other cosmological parameters. The result is shown as 
the first entry to Table [3] In the peak method we find the 
location of the peak position and use it as a standard ruler 
using Eq |15l This is a more direct, geometric test in the 
sense that the BAO is used strictly as a standard ruler (i.e. 
we are not using model-dependent shape of the correlation 
function), and in the sense that we constrain H{z) rather 
than integrals thereof. The results are shown as 3 bottom 
entries of Table O 

There are two parts to our error analysis. One is the er- 
rors on ^((T, n) itself, and the other is the errors on our mea- 
surement of the peak location. Both have been extensively 
tested with simulations. The error model for ^ is obtained 
from what is to our knowledge the largest ever cosmologi- 
cal simulation run to date, MICEL7680 with 453 Gpc^/h^, 
2048 ^ dark matter partic l es and 107 milhon h alos in a single 
box l|Fosalba et al.l l2008l : ICrocce et ai]|2009l ). We are able 
to create out of this simulation 216 independent mock LRG 
catalogs, each of which has a similar size to our SDSS DR6 
sample. Our mock catalogs have the same number density 
and very similar 2- and 3-point functions compared to data. 
Fig lAll shows, for the first time on BAO scales, a compari- 
son of the redshift space LOS linear model with non-linear 
measurements from simulations. Both model and simulation 
show very similar shapes and a prominent BAO peak which 
validates both our modeling and our simulated mocks. The 
error model constructed from these simulations is further 
validated by a comparison with jack-knife errors obtained 
from the data itself. Further details can be found in ij2.3l 
and Fig[5l 

The other part of our error analysis consists of using 
Monte Carlo simulations to simulate our peak measurement 
process. This is described in i^4.2l This method involves a 
further systematic uncertainty: because the data is quite 
noisy we can not be certain of what is the true model for 
the mean of the Monte Carlo simulations. We have used 
the data itself (or a smooth version of it) to measure both 
the statistical and systematic errors on these measurements. 
Thus, these errors are more model independent than the 
ones in the shape method, which depend strongly on the 
range of scales used for the fit. 

Given our measurements of H{z), it is natural to ask: 



H(z) rBAO + H„ HST + n„ LRG 




0.65 0,7 0.75 O.a 0.85 

o„,,= i-n^ 



Figure 22. The dotted line and continuous (black) open con- 
tours across the figure show the best fit and 1-sigma contour from 
q[z = 0.34] in Eq |18l ie using only our two independent H{z) 
measurements. Long-dashed line contour show Ax^ = 6.18 (2- 
sigma for 2 dof) constraints from combining our H(z) measure- 
ments with Ho from HST and our flm estimate from redshifts 
distortions (Paper I). Short-dashed contours (blue) show the cor- 
responding Ax^ = 6.18 constraints from the external data-set 
WMAP5, SNIa and BAO monopole. The inner red continuous 
contour shows Ax^ = 1 when combining both measurements. 

what constraint do they put on the dark energy abun- 
dance and equation of state? Our measurements of H{z) at 
z = 0.24 and z = 0.43 by themselves only weakly constrain 
f^DE and w. The constraint comes essentially from the accel- 
eration parameter q{z) in Ea llSI The mean and 1-sigma con- 
strains on S7de — w from q{z) are shown as dotted and con- 
tinuous black open contours across Fig |22l These constraints 
are very weak as compare d to the constra i nts fro m WMAP5 
-f SN -f monopole BAO iKomatsu et al.l (|2008l ). which are 
shown for Ax^ = 6.18 as a long-dashed (blue) line contour 
in Fig |22l The large short-dashed (black) line contour cor- 
responds to the 2-sigma constraints when we combine our 
H{z) measuremen ts with Hp = 70 ± 5 km/s/ Mpc, based on 
the HST analysis Ivan Leeuwen et "ahl l|2007l ). and our best 
fit fim = 0.245 ± 0.020 using the same LRG data in Paper 
I. The 2-sigma contours from the two independent data sets 
agree well, which is not trivial specially given the different 
assumptions and measurements involved in each data set. 
Combining both data sets, we find 

w = -0.957 ±0.053 (20) 
Qde = 0.734 ±0.023 (21) 

which is shown as red continuous contour in Fig |22l This 
is in excellent agreement with the cosmological constant 
model. The improvement in the error with respect to the 
1-sigma WMAP5±SN±BAO value oi w = -0.97 ± 0.06 
and is quite modest. A separate paper presents the impli- 
cation on cosmological parameters when using our radial 
BA O scale measurements in a way that is independent of 
Ho (jCaztanaga et al.ll2008l 'l. 

To illustrate what might be achievable in the near fu- 
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Figure 23. g(7r, cr) in data for the z=0.15-0.30 slice recalculated 
from real LRG positions with an added photo-z error of 0.003(1 + 
z) as predicted for the PAU survey. It should be compared with 
middle panel of Fig llOl which uses spectroscopic redshifts. 

ture, w e show in Fig. | 23 | t he expected ^ {a, n) for the PAU 
survey (|Benitez et al.ll2008l ) which has a photo-z accuracy of 
0.003(1 + z). We redo the analysis using current SDSS data 
(shce z=0.15-0.30) but we disperse the LRG positions by 
this photo-z error (we assume a Gaussian distribution with 
0.003(1-1-2) dispersion). The BAO signal does not appear to 
be washed away, as shown in Fig. 1231 This is expected be- 
cause the BAO peak is broader than the scale corresponding 
to the photo-z dilution. The PAU Survey proposes to map 
over 10 times the SDSS DR6 volume, ie to z=0.9 (this is 
a factor of 3 in sampling variance error) with 20 times the 
LRG number density (ie for L > L,) so shot noise will be 
negligible. All this should increase the signal-to-noise over 
SDSS DR6 by better than a factor of four. Thus, there is 
potential for a substantial improvement in the near future 
for the H{z) measurements using the techniques presented 
here. 
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Figure Al. Linc-of-sight (LOS) correlation averaged for < 
cr < 5 Mpc/h. Filled circles connected by lines represent the fol- 
lowing quantities measured from dark matter simulations: LOS 
mass autocorrelation in real space (blue dotted), redshift space 
(red dashed), cross real-redshift space (black solid). The dotted, 
dashed and solid lines (with no dots or circles) correspond to our 
analytic predictions for the LOS mass autocorrelation in real, red- 
shift and cross space. Finally, also shown is the monopole mass 
autocorrelation in real space (unfilled diamonds) which agrees 
very well with the prediction and, modulus sampling variance, 
should also agree with the LOS mass autocorrelation. 
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APPENDIX A: NUMERICAL 

INVESTIGATIONS OF NONLINEAR EFFECTS 

In this Appendix, we describe a number of measurements we 
have made to gauge the size of various nonlinear effects on 
the magnification bias correction to the galaxy correlation. 
The measurements are done on numerical N-body simula- 
tions described in the main text. Recall that the magnifica- 
tion bias correction to galaxy overdensity at radial distance 
X2 is given by line of sight integral: 



K 



dx 



(X2 - x)x 

a(x)X2 



5{X) 



(Al) 



where K ^ {^s - 2)§HiQ„ 



Gonsider now a regrouping of 



the terms in Eq.(l7l: Sots = Sg + 5 fj,, where 5g corresponds 
now to the galaxy overdensity in redshift space (i.e. includ- 
ing non- linear redshift space distortions). The cross-term 
contribution to the observed correlation function is then 
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Figure A2. Magnification bias correction based on Eq. I IA2I I 
(which replaces Eq. [9]) for a mean redshift z = 0.34, Qrn = 0.25 
and s = 1.5 for Hnear (bottom lines) or non-linear (upper lines) 
bias using: a) the prediction of the cross-correlation ggm in real 
space (dotted), b) the prediction of the cross-correlation ^gm in 
cross space (solid). Here, cross space means the galaxy is in red- 
shift space and the mass is in real space. This Jg^i is computed 
with the transverse separation cr averaged over — 5 Mpc/h. 



= {Sg{l)5^{2)) , where 1 and 2 denotes the positions 
of the the foreground and background galaxies. The LOS 
separation in redshift space is tt, and the transverse sepa- 
ration is a. The relevant galaxy-magnification correlation is 
(replacing Eq. [5]): 



K 



, (X2 - x)x g ( 



x,<y) (A2) 



where = {5aiX)5m{x)) is the cross correlation of galaxy 
fluctuations in redshift space with dark matter fluctuations 
in real space. Here, Xi is the redshift space radial distance to 
galaxy 1, but ^2 and x are real space distances. Note that 
^gm depends only on the pair separation i.e. (xi — XiO"); 
while ^gp depends both on X2 and on the pair separation 
(7r,a). 

In our new calculation here we will fix X2 to be at the 
mean redshift of the survey, but we have also studied how the 
result changes when we take into account the fact that X2 
varies across the survey (see below) - basically, if xi are X2 
are well separated, ^g,t is almost independent of X2 (except 
indirectly through the LOS separation tt). One can see how 
Eq. ((9| comes about from Eq. (|A2|I - take the limit in which 
X2 — Xi ^ X2,Xii and ignore the subtle difference between 
redshift space and real space distances. Our task here is to 
evaluate Eq. (|A2[I exactly using simulations. 

First, as a warm-up, we use the dark matter (as op- 
posed to halo) simulations which have the correct veloc- 
ities on small scales. We compare the mass-mass correla- 
tion ^mm(7r,(T) iu real space, redshift space and the cross- 
correlation between real and redshift space. (Note: S^mm can 
be thought of as a first proxy to £,gm\ we will examine ^haio,m 
below as a better proxy.) We smooth ^mm(7r,o") in bins of 
5Mpc/h, as we do with real LRG galaxies. In Fig lAll we show 



a set of N-body measurements (shown as dots connected by 
lines) for the dark matter LOS auto-correlation in real space 
(dotted blue line), redshift space (dashed red line) and in 
the cross-correlation of real and redshift space (solid black 
line). Note how the BAO peak in redshift or cross-space is 
enhanced quite a bit compared to that in real space. We 
compare these simulation results with analytic predictions 
(lines without dots; blue dotted for real space, red dashed 
for redshift space and black solid for cross-space) which com- 
bine linear theory on large scales with velocity dispersion 
on small scales (see paper I for an extensive explanation of 
the model). Note how the analytic model roughly matches 
the simulation data. However, despite the similarities in the 
shape, the amplitude of the LOS baryon peak seems slightly 
larger in the simulations than in the analytic model. It is 
not entirely clear to us how significant the difference is, but 
it should be noted that non-linear redshift space distortions 
have yet to be proper l y exp lored and modeled in the LOS 
direction IScoccimarrol ( 2004 ) . We will see a similar tendency 
in the actual LRG galaxies but the measurement errors are 
large. Ultimately, what we need is something like a LOS 
integral of (see Eq. |A2j ). We find that the difference 
between integrating 5mm(7i",(j) in real space, redshift space 
or in the cross space is not very significant. 

Next, we repeat the same exercise for ^gm, using halos 
in the simulations as a proxy for LRGs. The conclusions are 
quite similar to those for ^mm. In particular, we find that 
as far as the LOS integral of is concerned (Eq. [A2) ). it 



actually does not matter much whether we use ^gm in real 
space, redshift space or cross-space. What distinguishes ^gm 
from ^mm is the galaxy bias which is about 2 on large scales 
but rise to about 3 — 4 on small scales. As emphasized in 
H2.3I the nonlinear galaxy bias affects C9m('''"i f^'' a small 
cr, even if tt is large. This is illustrated in Fig. IA2I The 
correct ^g^ using a nonlinear galaxy bias is larger than the 
one using a linear galaxy bias by about a factor of 1.5. Note 
also how the n dependence is close to linear, supporting the 
approximations made in Eq. 

Finally, we have also studied the effect of the chosen 
value for X2 in the integral of Ea lA2l The real magnification 
should be a weighting of integrals at the different redshifts 
(or X2) of the survey. Each redshift bin is weighted depend- 
ing on the number of galaxies and the volume under this X2. 
The exact result is almost identical to taking X2 to be at the 
median redshift of the survey. 
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